Local transport in a disorder-stabilized correlated insulating phase 
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We report the experimental realization of a correlated insulating phase in 2D GaAs/AlGaAs 
heterostructures at low electron densities in a limited window of background disorder. This has 
been achieved at mesoscopic length scales, where the insulating phase is characterized by a universal 
hopping transport mechanism. Transport in this regime is determined only by the average electron 
separation, independent of the topology of background disorder. We have discussed this observation 
in terms of a pinned electron solid ground state, stabilized by mutual interplay of disorder and 
Coulomb interaction. 

PACS numbers: 73.21.-b, 73.20.Qt 



In the presence of Coulomb interaction, both mag- 
netic field and disorder are predicted to stabilize many- 
body charge-ordered ground states Strong perpendic- 
ular magnetic field quenches the vibrational motion 
of electrons, and has been extensively exploited to real- 
ize a charge-density wave (CDW) ground state in sys- 
tems with weak background disordcr.2i4 Despite the ef- 
fort however, the nature of localization in such systems 
has been controversial, with both pinned Wigner solid 
(WS) formation and inhomogeneity-driven percolation 
transition being suggested^ On the other hand, disor- 
der stabilizes Coulomb correlation effects by introducing 
a pinning gap Apin in the phonon density of states, which 
provides a long wavelength cutoffi^ This has led to the 
theoretical prediction of several forms of CDW ground 
states at zero or low B±^ . Systematic experimental inves- 
tigations on such possibilities, however, have been rare, 
and form the subject of this work. 

Increasing the magnitude of background potential 
fluctuations increases Apin which stabilizes the CDW 
phases to higher temperatures. In modulation-doped 
GaAs/AlGaAs heterostructures, where disorder primar- 
ily arises from the charged dopant ionS)& Apin ^ 
exp (— 47r(5sp/ -s/S^ee) depends strongly on the setback dis- 
tance (5sp that separates the 2D electron system (2DES) 
and the dopants, where is the mean distance between 
the electrons in the 2DESi^ However, disorder affects 
the ground state transport in two critical ways. First, 
presence of Apin disintegrates the CDW phase into do- 
mains of finite size Ad ^ sound velocity/Apin. At strong 
pinning. Ad becomes microscopically small, leading to sig- 
nificant averaging in transport measurements with con- 
ventional macroscopic devices. Second, strong potential 
fluctuations can also result in a "freezing" of transport 
below a certain percolation threshold even when electron 
density (rig) is relatively high, thereby making the regime 
of strong effective Coulomb interaction inaccessible. 

Here, we show that these difficulties can be largely 
overcome by using modulation-doped heterostructures of 
mesoscopic dimensions. In such devices transport freezes 
at much lower rig in comparison to macroscopic devices 
at the same 5sp or disorder, thereby allowing transport at 
a large interaction parameter rg = l/ag^Trrig '--^ 7— 8 (og 



is the effective Bohr radius), even when 5s-p is relatively 
small. Typical dimension L of our devices in the current 
carrying direction was chosen to be ~ 2 — 4^m, which 
is also similar in order of magnitude to the Ad suggested 
by recent microwave absorption studies for pinned WS 
ground states^ The low-i?j^ magnetotransport in these 
devices was found to display a striking universality in 
that the hopping distance in the localized regime was de- 
termined by = 1/ ^/n^^ rather than the details of back- 
ground disorder, indicating an unusual self-localization of 
electrons at sufficiently low rig. 

We have used Si modulation-doped GaAs/AlGaAs het- 
erostructures where (5gp was varied from 20 nm to 80 
nm. At a fixed ng, the effect of (5gp on the strength of 
potential fluctuations is reflected in the mobility /i, as 
can be observed from Fig. lb. Both monolayer [5)- and 
bulk-doped wafers were used. Relevant properties of the 
devices are given in Table I. Devices were cooled from 
room temperature to 4.2 K over 24 - 36 hours to allow 
maximal correlation in the dopant layer (redistribution 
of DX-centers)iS This slow cooldown technique also lead 
to excellent reproducibility over repeated thermal cycles. 
Electrical measurements were carried out with standard 
low frequency (7.2 Hz) four-probe technique with exci- 
tation current of ~ 0.01 — 0.1 nA to minimize heating 



TABLE I: Geometrical and structural property of the de- 
vices, ns is the density of Si dopants, and W is the width. 
The background doping concentration is < lO^^cm^'' in all 



devices. 
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"The doping concentration of bulk doped devices is 2 X lO^^cm 
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FIG. 2: Typical magnetoresistivity traces in four samples 
with varying level of disorder. The vertical lines denote 
v — 1. The numbers indicate electron density in units of 
10^*^ cm^^. Be denotes the field scale up to which a quadratic 
Bx-dependence could be observed. The parameters a and pe 
were obtained from the slope and y-intercept of linear fits to 
ln(p) - Bitr aces, respectively. 



FIG. 1: (Color online) (a) Conductivity (a) of mesoscopic 
samples as a function of electron density ns at T « 0.3 K. The 
crosses denote n* for individual samples (see text). Inset: rig- 
dependence of a for a macroscopic Hall bar A77L. The solid 
line is the best fit of a classical percolation-like scaling relation 
a ^ (ris — ric)^. (b) Ssp dependence of mobility at constant Us 
and ns for heterostructures similar to those used in presented 
work (c) Resistivity (p) as a function of temperature measured 
at B± = 1 T. The solid line represent a power law of ~ T~^, 
vertical lines in (c) and (d) indicate the Fermi temperatures 
Tp. (d) Activation and saturation of p at B± — 1.5 T. 



and other nonlinear effects. A direct measurement of Ug 
within the microscopic region was carried out with an 
edge-state reflection-based techniquciiS 

In Fig. la we compare the ng-scale of localization tran- 
sition at B± = and T = 0.3 K in macroscopic and 
microscopic devices from the same wafer. In a standard 
100/.tmx900/im Hall bar, as illustrated with wafer A2677, 
the linear conductivity a — > (A77L: Inset of Fig. la) 
at ~ 3 times the rig compared to the mesoscopic sample 
(A77) from the same wafer. Further, a in the large sam- 
ple A77L shows excellent classical percolation-like scaling 
a ~ (rig — Tic)''' (ric = 1.72 x 10^" cm^^), where 7 « 2, im- 
plying a inhomogeneity-driven percolation transition at 
nonzero T— (solid line in inset of Fig. la). Similar scaling 
in the mesoscopic systems, however, was found to be dif- 
ficult, with unphysically large estimates of 7 ^ 3.2 — 3.7 
(not shown), indicating a different mechanism of local- 
ization transition. 

As rig is lowered below a sample-dependent character- 
istic scale n* (denoted by the crosses in Fig. la), onset 
of strong localization is identified by the resistivity p (= 
1/cr) exceeding « h/e'^. At Ug n*, the T-dependence 



of p at a fixed rig can be divided into three regimes, as 
illustrated with A78a: First, transport in the classical 
regime at T > Tp is magnified in Fig. Ic, where Tp is 
the Fermi temperature. In this regime p oc , where 
/3 ~ 1 (indicated by the sohd line). As T is decreased, the 
onset of the quantum regime (T < Tp) results in stronger 
increase in p with decreasing T. Note that the clear clas- 
sical to quantum crossover implies a well-defined Tp, and 
hence a uniform charge density distribution down to the 
lowest Hs ~ 6.5 X 10® cm~^ (In A77L, inhomogcneity sets 
in at Us as large as ^ 4 — 5 x 10^° cm^^). In the quan- 
tum regime and for > T > T* , Fig. Id shows that the 
behavior of p is activated with p{T) = pa exp(e3/fcBr), 
where £3 is the activation energy. From the Ug- and B±- 
depcndence of the pre-exponential pa , we have shown ear- 
lier that the transport mechanism in this regime corre- 
sponds to nearest-neighbor hopping>i£ Below the char- 
acteristic scale T* ~ 1 K, variation of p becomes weak, 
tending to a finite magnitude even in the strongly lo- 
calized regime. This saturation in the insulating regime 
cannot be explained in terms of an elevated electron tem- 
perature due to insufficient thermal coupling to the lat- 
tice since T* depends only weakly on electron density up 
to Us ^ n* (Fig. Id), and the damping of Shubnikov-de 
Haas oscillations in the metallic regime shows the base 
electron temperature to be « 300 mK. 

In order to explore the physical mechanism behind the 
weak T-dependence of p, we have carried out extensive 
magnetoresistivity (MR) measurements at the base T. 
Figs. 2a to 2d show the B^-dependence of MR in the 
insulating regime of four devices with increasing Ssp from 
20 nm to 80 nm. In general, we find a strong negative MR 
in A07, A78 and C67 at low B±, which can be attributed 
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FIG. 3: Absolute magnitude of a obtained from the slope of 
ln(p) — B\ traces for five different samples. Inset shows the 
same data in a log-log scale. The slope of the solid line is 
-3/2. 



to interference of hopping paths. The negative MR is 
followed by an exponential rise in p as B±^ is increased 
further. We have recently shown that the logarithm of 
such a positive MR at low B± varies in a quadratic man- 
ner with B^, i.e., p{B±) = pB exp(ai3^), where pB and a 
are rig-dependent factors liS Such a variation, denoted by 
the solid lines in Fig. 2, is found to be limited to rig '^n*, 
and extends over a i3j^-scalc of Be, where Be was found 
to decrease rapidly as ^sp is increased. Note that in T46 
(lowest disorder), neither a clear negative MR nor an 
exponential B^-dependence were observed. A physical 
significance of Be and of the qualitatively different MR 
behavior of T46 will be discussed later. 

The observed behavior of p can be naturally explained 
in the framework of tunnelling of electrons between two 
trap sites separated by a distance ry . In weak B± , such 
that the magnetic length A = \/h/eB_i_ ^ ^, where ^ is 
the localization length, the asymptotic form of the hy- 
drogenic wave function changes from ip{r) ^ exp(— r/^) 
to tpir) - exp(-r/^ - r3^/24A^)pii This leads to a MR, 
p(Bj_) = pQCxp{2nj/C)cxp{Ce^rf^^Bl/12h'^), which im- 
plies 



PB = Po exp(2rij / ^) and 



ie/12ft^ (1) 



While pb depends on the tunnelling probability at B± ~ 
0; a denotes the rate of change of this probability when 
B± is switched on. Importantly, both parameters provide 
information on the intersitc distance , as well as ^ in- 
dependently. The parameter C ^ 0.5 — 1 depends on the 
number of bonds at percolation threshold in the random 
resistor network (we shall subsequently assume C ~ 1). 
Since conventional hopping sites are essentially impurity 
states, both a and pB are expected to be strongly dis- 
order dependent. Note that, since wave function overlap 
plays a critical role in transport, a direct source-to-drain 
tunnelling is ruled out in our case^ 



From the MR data we have evaluated a and pB from 
the slope and intercept of the ln{p) — B\ traces. Fur- 
ther details of the analysis can be found elsewherepi^i In 
Fig. 3 we have shown a as a function of ng for five dif- 
ferent samples up to the corresponding n*. Strikingly, 
the absolute magnitudes of a from different samples are 
strongly correlated, and can be described by an universal 
ng-depcndent function over nearly two orders of magni- 
tude. At stronger disorder (e.g., A07), locahzation occurs 
at a higher rig resulting in a lower a, while at lower dis- 
order (e.g. C67) localization occurs at lower rig yielding 
a larger magnitude of a. This indicates that magneto- 
transport in such mesoscopic samples is not determined 
directly by disorder, but by rig in the localized regime. 
Qualitatively, the decreasing behavior of a with increas- 
ing rig itself is inconsistent with the single-particle local- 
ization in an Anderson insulator^S*i^ 

From the strong sample-to-sample correlation in the 
magnitude of a, a disorder-associated origin of nj , is 
clearly unlikely. For example, taking r-^ ^ Jgp will lead 
to distinct sets of a for wafers with different Ss-p- How- 
ever, in the context of a pinned CDW ground state, 
another relevant length scale is reo. Indeed, in case of 
tunnelling events over a mean electron separation, i.e., 
Hj ~ ?'cc, we find that Eq.ndescribes both absolute mag- 
nitude, as well as the rig-dependence of a quantitatively. 
Using Tij « Eq. pleads to a cx rig as in- 

deed observed experimentally (solid line in the inset of 
Fig. 3). Allowing for sample-to-sample variation, we find 

3/2 
Is 

we get ^ = 9.0 ± 2.6 nm 
(« 10.5 nm). 

The analysis can be immediately checked for consis- 
tency from the roo-dependence of pb- From Fig. 4, we 
find that pe increases strongly with increasing r^c when 
rig <C n*, as expected in the simple tunnelling frame- 
work (Eq. In spite of the scatter, the overall slopes 
of the In(pB) — ''ee plots are similar in different samples 
(solid lines) with S, estimated to be w 13 ± 4 nm, agree- 
ing with that obtained from the analysis of a. Note that 
the Pb deviates from the exponential dependence, and 
tends to saturate as rig n*. While this is not com- 
pletely understood at present, we note that the satura- 
tion in Pb occurs within the range pe ^ 1 ~ 2 x /i/e^, 
irrespective of sample details. Similar universality in the 
hopping pre-exponcntial has been observed in the con- 
text of T-dependence of p in variable-range hopping^ 
and has been suggested to indicate an electron-electron 
interaction mediated energy-transfer mechanism. 

We now discuss the physical scenario which could lead 
to the electron separation-dependent hopping transport. 
We show that our observations can be explained in the 
theoretical framework of defect motion in a quantum 
solid that was originally developed by Andreev and Lif- 
shitz in context of solid He^ ^ and later adapted for a WS 
ground state^Siii In our case, transport in both quan- 
tum and classical regime can be understood in terms of 
tunnelling of localized defects in an interaction-induced 



a = (1.7 ±0.5) X lO^i/rig^^ T'^ from which, using Eq.Cl 

which is close to in GaAs 
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FIG. 4: (Color online) The dependence of ps on the average 
electron separation ree in five different samples. The slope of 
the solid lines gives an estimate of ^ (Eq.Ql. 

pinned electron solid phase as ng is reduced below the 
melting point n*. The defects, which act as quasipar- 
ticles at low T, can arise from regular interstitials, va- 
cancies, dislocation loops etc., as well as from zero-point 
vibration of individual lattice sitesii^ The scale of zero- 
point fluctuation ~ h/ree^/m*Uc ~ ^ir/^yr^ > 1, is in- 
deed strong in our case over the experimental range of 
rig, where J7c ~ e^/47rerco is the interatomic interaction 
energy scale. 

In the quantum regime, the transport at higher T 
(Tp » T 3> r*) is predicted to be thermally activated 
nearest- neighbor hopping of localized defects, while at 
lower T (<C T*) tunnelling of such defects leads to a T- 
independent transport pi^ While this clearly describes the 
weak T-dependence of p at low temperatures (Fig. Id), 
the strongest support to this picture comes from the fact 
that the natural length scale of tunnelling is indeed the 
average electron separation roc. This immediately ex- 
plains the unusual Ug (or roe)-dependence of both a and 
PBj as well as the apparent insensitivity of these parame- 
ters to local disorder. The negative MR at low B± caused 
by destruction of interference is then expected to persist 



up to a B± corresponding to v = nsh/eB± ~ 1 (one 
flux quantum 0o within an area of r^g), as indeed ob- 
served in our experiments (Fig. 2). The tunnelling of 
defect scenario also allows an estimate of the crossover 
scale kBT* = £3/ ln(Apin/Ae)fi^ where Ae is the band- 
width. For a pinned WS ground state, using the expres- 
sion of Apin in Ref. [8], experimentally measured £3, and 
Ae ft,^/8TO*r^(,, we find T* ~ 0[1 Kelvin] over the 
experimental range of Ug in A78a, giving good order-of- 
magnitude agreement to the observed scale of T*. Fi- 
nally, the behavior of p ~ T^^ in the classical regime 
(T > TF)(Fig. Ic) has also been recently observedji^ and 
interpreted in terms of transport mediated by defect- 
type topological objects (Fermi-liquid droplets) in the 
WS phaseia 

In presence of pinning, the MR-data suggests the 
asymptotic form of the wave function 'ijj{r) ^ cxp(r/^), 
where ^ ~ Og. However, the interplay of confinement 
arising from the magnetic potential and disorder pinning 
is expected to be critical in determining ■0(r), with dis- 
order pinning dominating at low B± . This is expected to 
result in the upper cutoff Be that decreases with decreas- 
ing disorder, as observed experimentally. The intricate 
interplay between disorder, electron-electron interaction 
and magnetic field is further illustrated by the absence 
of a clear B^ dependence of the MR in T46 (largest Ssp), 
which could be explained by a prohibitively small B^ or 
the very instability of the solid phase at sufhcicntly low 
disorder. On the other hand, devices with 5sp < lOnm 
showed inhomogeneity-driven Coulomb blockade oscilla- 
tions in the localized regime, making the investigation of 
such a charge correlated state impossible. A quantitative 
understanding of the scale of Be, as well as the specific 
spatial structure of the ground state in the intermediate 
disorder regime, will require further investigations, which 
are presently in progress. 
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